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ABSTRACT 

We present new full- sky temperature and polarization maps in five frequency bands from 23 
to 94 GHz, based on data from the first five years of the WMAP sky survey. The new maps are 
consistent with previous maps and are more sensitive. The five-year maps incorporate several 
improvements in data processing made possible by the additional years of data and by a more 
complete analysis of the instrument calibration and in-flight beam response. We present several 
new tests for systematic errors in the polarization data and conclude that W band polarization 
data is not yet suitable for cosmological studies, but we suggest directions for further study. We 
do find that Ka band data is suitable for use; in conjunction with the additional years of data, 
the addition of Ka band to the previously used Q and V band channels significantly reduces the 
uncertainty in the optical depth parameter, r. Further scientific results from the five year data 
analysis are presented in six companion papers and are summarized in §7 of this paper. 

With the 5 year WMAP data, we detect no convincing deviations from the minimal 6- 
parameter ACDM model: a flat universe dominated by a cosmological constant, with adiabatic 
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and nearly scale-invariant Gaussian fluctuations. Using WMAP data combined with measure- 
ments of Type la supernovae (SN) and Baryon Acoustic Oscillations (BAO) in the galaxy dis- 
tribution, we find (68% CL uncertainties): Sl b h 2 — 0.02265 ± 0.00059, f2 c /i 2 = 0.1143 ± 0.0034, 
n A = 0.721 ± 0.015, n s = 0.9601q.oi3, t = 0.084 ± 0.016, and A|, - (2.4571^1) x HT 9 at 
k = 0.002 Mpc -1 . From these we derive: <7g = 0.817 ± 0.026, Ho = 70.1 ± 1.3 km s -1 Mpc -1 , 
n b = 0.0462 ± 0.0015, n c = 0.233 ± 0.013, = 0.1369 ± 0.0037, z rei0 n = 10.8 ± 1.4, and 

to = 13.73 ±0.12 Gyr. The new limit on the tensor-to-scalar ratio is r < 0.20 (95% CL), while the 
evidence for a running spectral index is insignificant, dn s /dlnk = — 0. 032^ 020 (68% CL). We 
obtain tight, simultaneous limits on the (constant) dark energy equation of state and the spatial 
curvature of the universe: —0.11 < 1+tu < 0.14 (95% CL) and —0.0175 <Qk < 0.0085 (95% CL). 
The number of relativistic degrees of freedom, expressed in units of the effective number of neu- 
trino species, is found to be jV e ff = 4.4 ±1.5 (68% CL), consistent with the standard value of 
3.04. Models with N e ft = 0 are disfavored at >99.5% confidence. Finally, new limits on phys- 
ically motivated primordial non-Gaussianity parameters are —9 < /]y£ al < 111 (95% CL) and 
— 151 < f e *™ 1 < 253 (95% CL) for the local and equilateral models, respectively. 

Subject headings: cosmic microwave background, cosmology: observations, early universe, dark 
matter, space vehicles, space vehicles: instruments, instrumentation: detectors, telescopes 


1. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) is a Medium-Class Explorer (MIDEX) satellite 
aimed at elucidating cosmology through full-sky observations of the cosmic microwave background (CMB). 
The WMAP full-sky maps of the temperature and polarization anisotropy in five frequency bands provide 
our most accurate view to date of conditions in the early universe. The multi-frequency data facilitate the 
separation of the CMB signal from foreground emission arising both from our Galaxy and from extragalactic 
sources. The CMB angular power spectrum derived from these maps exhibits a highly coherent acoustic 
peak structure which makes it possible to extract a wealth of information about the composition and history 
of the universe, as well as the processes that seeded the fluctuations. 

WMAP data have played a key role in establishing ACDM as the new standard model of cosmology 
(Bennett et al. 2003; Spergel et al. 2003; Hinshaw et al. 2007; Spergel et al. 2007): a flat universe dominated 
by dark energy, supplemented by dark matter and atoms with density fluctuations seeded by a Gaussian, 
adiabatic, nearly scale invariant process. The basic properties of this universe are determined by five numbers: 
the density of matter, the density of atoms, the age of the universe (or equivalently, the Hubble constant 
today), the amplitude of the initial fluctuations, and their scale dependence. 

By accurately measuring the first few peaks in the angular power spectrum, WMAP data have enabled 
the following accomplishments: 

• Showing the dark matter must be non-baryonic and interact only weakly with atoms and radiation. 
The WMAP measurement of the dark matter density puts important constraints on supersymmetric 
dark matter models and on the properties of other dark matter candidates. With five years of data 
and a better determination of our beam response, this measurement has been significantly improved. 

• Precise determination of the density of atoms in the universe. The agreement between the atomic 
density derived from WMAP and the density inferred from the deuterium abundance is an important 
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test of the standard big bang model. 

• Determination of the acoustic scale at redshift z — 1090. Similarly, the recent measurement of baryon 
acoustic oscillations (BAO) in the galaxy power spectrum (Eisenstein et al. 2005) has determined 
the acoustic scale at redshift z ~ 0.35. When combined, these standard rulers accurately measure 
the geometry of the universe and the properties of the dark energy. These data require a nearly flat 
universe dominated by dark energy consistent with a cosmological constant. 

• Precise determination of the Hubble Constant, in conjunction with BAO observations. Even when 
allowing curvature (Qq ^ 1) and a free dark energy equation of state (w ^ —1), the acoustic data 
determine the Hubble constant to within 3%. The measured value is in excellent agreement with inde- 
pendent results from the Hubble Key Project (Freedman et al. 2001), providing yet another important 
consistency test for the standard model. 

• Significant constraint of the basic properties of the primordial fluctuations. The anti-correlation seen 
in the temperature/polarization (TE) correlation spectrum on 4° scales implies that the fluctuations 
are primarily adiabatic and rule out defect models and isocurvature models as the primary source of 
fluctuations (Peiris et al. 2003). 

Further, the WMAP inference of the primordial power spectrum of matter fluctuations constrains the 
physics of inflation, our best model for the origin of these fluctuations. Specifically, WMAP data provide 
the best measurement to date of the scalar spectrum’s amplitude and slope, and place the most stringent 
limits to date on the amplitude of tensor fluctuations. 

The statistical properties of the CMB fluctuations measured by WMAP are close to Gaussian; however, 
there are several hints of possible deviations from Gaussianity. Significant deviations would be a very 
important signature of new physics in the early universe. 

Large-angular-scale polarization measurements currently provide our best window into the universe at 
z ~ 10. The WMAP data imply that the universe was reionized long before the epoch of the oldest known 
quasars. By accurately constraining the optical depth of the universe, WMAP not only constrains the age 
of the first stars but also determines the amplitude of primordial fluctuations to better than 3%. This result 
is important for constraining the growth rate of structure. 

This paper summarizes results compiled from 5 years of WMAP data that are fully presented in a suite 
of 7 papers (including this one). The new results improve upon previous results in many ways: additional 
data reduces the random noise, which is especially important for studying the temperature signal on small 
angular scales and the polarization signal on large angular scales; five independent years of data enable 
comparisons and null tests that were not previously possible; the instrument calibration and beam response 
have been much better characterized, due in part to improved analyses and to additional years of data; and, 
other cosmological data have become available. 

In addition to summarizing the other papers, this paper reports on changes in the WMAP data process- 
ing pipeline, presents the 5 year temperature and polarization maps, and gives new results on instrument 
calibration and on potential systematic errors in the polarization data. Hill et al. (2008) discuss the program 
to derive an improved physical optics model of the WMAP telescope, and use the results to better determine 
the WMAP beam response. Gold et al. (2008) present a new analysis of diffuse foreground emission in the 
WMAP data and update previous analyses using 5 year data. Wright et al. (2008) analyze extragalactic 
point sources and provide an updated source catalog, with new results on source variability. Nolta et al. 
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(2008) derive the angular power spectra from the maps, including the TT, TE, TB, EE, EB, and BB spectra. 
Dunkley et al. (2008) produce an updated likelihood function and present cosmological parameter results 
based on 5 year WMAP data. They also develop an independent analysis of polarized foregrounds and use 
those results to test the reliability of the optical depth inference to foreground removal errors. Komatsu et al. 
(2008) infer cosmological parameters by combining 5 year WMAP data with a host of other cosmological 
data and discuss the implications of the results. Concurrent with the submission of these papers, all 5 year 
WMAP data are made available to the research community via NASA’s Legacy Archive for Microwave Back- 
ground Data Analysis (LAMBDA). The data products are described in detail in the WMAP Explanatory 
Supplement (Limon et al. 2008), also available on LAMBDA. 

The WMAP instrument is composed of 10 differencing assemblies (DAs) spanning 5 frequencies from 23 
to 94 GHz (Bennett et al. 2003): 1 DA each at 23 GHz (Kl) and 33 GHz (Kal), 2 each at 41 GHz (Q1,Q2) 
and 61 GHz (V1,V2), and 4 at 94 GHz (W1-W4). Each DA is formed from two differential radiometers 
which are sensitive to orthogonal linear polarization modes; the radiometers are designated 1 or 2 (e.g., VI 1 
or W12) depending on polarization mode. 

In this paper we follow the notation convention that flux density is S ~ and antenna temperature 
is T ~ where the spectral indices are related by /3 = a — 2. In general, the CMB is expressed in 
terms of thermodynamic temperature, while Galactic and extragalactic foregrounds are expressed in antenna 
temperature. Thermodynamic temperature differences are given by AT = A T^[(e x - l) 2 /;r 2 e x ], where 
x = hv/kTo , h is the Planck constant, v is the frequency, k is the Boltzmann constant, and T 0 = 2.725 K is 
the CMB temperature (Mather et al. 1999). A WMAP band-by-band tabulation of the conversion factors 
between thermodynamic and antenna temperature is given in Table 1. 


2. CHANGES IN THE 5 YEAR DATA ANALYSIS 

The 1 year and 3 year data analyses were described in detail in previous papers. In large part, the 5 year 
analysis employs the same methods, so we do not repeat a detailed processing description here. However, we 
have made several improvements that are summarized here and described in more detail later in this paper 
and in a series of companion papers, as noted. We list the changes in the order they appear in the processing 
pipeline: 

• There is a ~ 1' temperature-dependent pointing offset between the star tracker coordinate system 
(which defines spacecraft coordinates) and the instrument boresights. In the 3 year analysis we in- 
troduced a correction to account for the elevation change of the instrument boresights in spacecraft 
coordinates. With additional years of data, we have been able to refine our thermal model of the 
pointing offset, so we now include a small (<1') correction to account for the azimuth change of the 
instrument boresights. Details of the new correction are given in the 5 year Explanatory Supplement 
(Limon et al. 2008). 

• We have critically re-examined the relative and absolute intensity calibration procedures, paying special 
attention to the absolute gain recovery obtainable from the modulation of the CMB dipole due to 
WMAP’ s motion. We describe the revised procedure in §4 and note that the sky map calibration 
uncertainty has decreased from 0.5% to 0.2%. 

• The WMAP beam response has now been measured in 10 independent “seasons” of Jupiter observa- 
tions. In the highest resolution W band channels, these measurements now probe the beam response 



~44 dB down from the beam peak. However, there is still non-negligible beam solid angle below this 
level (~0.5%) that needs to be measured to enable accurate cosmological inference. In the 3 year 
analysis we produced a physical optics model of the A-side beam response starting with a pre-flight 
model and fitting in-flight mirror distortions to the flight Jupiter data. In the 5 year analysis we 
have extended the model to the B-side optics and, for both sides, we have extended the fit to include 
distortion modes a factor of 2 smaller in linear scale (4 times as many modes). The model is used to 
augment the flight beam maps below a given threshold. The details of this work are given in Hill et al. 
(2008). 

The far sidelobe response of the beam was determined from a combination of ground measurements 
and in-flight lunar data taken early in the mission (Barnes et al. 2003). For the current analysis, we 
have replaced a small fraction of the far sidelobe data with the physical optics model described above. 
We have also made the following changes in our handling of the far sidelobe pickup (Hill et al. 2008): 
1) We have enlarged the “transition radius” that defines the boundary between the main beam and 
the far sidelobe response. This places a larger fraction of the total beam solid angle in the main beam 
where uncertainties are easier to quantify and propagate into the angular power spectra. 2) We have 
moved the far sidelobe deconvolution into the combined calibration and sky map solver (§4). This 
produces a self-consistent estimate of the intensity calibration and the deconvolved sky map. The 
calibrated time-ordered data archive has had an estimate of the far sidelobe response subtracted from 
each datum (as it had in the 3 year processing). 

We have updated the optimal filters used in the final step of map-making. The functional form of the 
filter is unchanged (Jarosik et al. 2007), but the fits have been updated to cover years 4 and 5 of the 
flight data. 

Each WMAP differencing assembly consists of two radiometers that are sensitive to orthogonal linear 
polarization states. The sum and difference of the two radiometer channels split the signal into intensity 
and polarization components, respectively. However, the noise levels in the two radiometers are not 
equal, in general, so more optimal sky map estimation is possible in theory, at the cost of mixing 
intensity and polarization components in the process. For the current analysis, we investigated one such 
weighted algorithm and found that the polarization maps were subject to unacceptable contamination 
by the intensity signal in cases where the beam response was non-circular and the gradient of the 
intensity signal was large, e.g., in K band. As a result, we reverted to the unweighted (and unbiased) 
estimator used in previous work. 

We have improved the sky masks used to reject foreground contamination. In previous work, we defined 
masks based on contours of the K band data. In the 5 year analysis we produce masks based jointly 
on K band and Q band contours. For a given sky cut fraction, the new masks exclude flat spectrum 
(e.g. free-free) emission more effectively. The new masks are described in detail in Gold et al. (2008) 
and are provided with the 5 year data release. In addition, we have modified the “processing” mask 
used to exclude very bright sources during sky map estimation. The new mask is defined in terms 
of low-resolution (r4) sky pixels to facilitate a cleaner definition of the pixel-pixel inverse covariance 
matrices, AT -1 . One side effect of this change is to introduce a few r4-sized holes around the brightest 
radio sources in the analysis mask, which incorporates the processing mask as a subset. 

We have amended our foreground analysis in the following ways: 1) Gold et al. (2008) perform a pixel- 
by-pixel analysis of the joint temperature and polarization data to study the breakdown of the Galactic 
emission into physical components. 2) We have updated some aspects of the Maximum Entropy (MEM) 
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based analysis, as described in Gold et al. (2008). 3) Dunkley et al. (2008) develop a new analysis of 
polarized foreground emission using a Gibbs sampling approach that yields a cleaned CMB polarization 
map and an associated covariance matrix. 4) Wright et al. (2008) update the WMAP point source 
catalog and present some results on variable sources in the 5 year data. However, the basic cosmological 
results are still based on maps that were cleaned with the same template-based procedure that was 
used in the 3 year analysis. 

• We have improved the final temperature power spectrum, Cj T ', by using a Gibbs-based maximum 
likelihood estimate for l < 32 (Dunkley et al. 2008) and a pseudo -Ci estimate for higher / (Nolta et al. 
2008). As with the 3 year analysis, the pseudo-C/ estimate uses only V- and W-band data. With 5 
individual years of data and six V- and W-band differencing assemblies, we can now form individual 
cross-power spectra from 15 DA pairs within each of 5 years and from 36 DA pairs across 10 year pairs, 
for a total of 435 independent cross-power spectra. 

• In the 3 year analysis we developed a pseudo-C* method for evaluating polarization power spectra in 
the presence of correlated noise. In the present analysis we additionally estimate the TE, TB, EE, EB, 
& BB spectra and their errors using an extension of the maximum likelihood method in Page et al. 
(2007). However, as in the 3 year analysis, the likelihood of a given model is still evaluated directly 
from the polarization maps using a pixel-based likelihood. 

• We have improved the form of the likelihood function used to infer cosmological parameters from the 
Monte Carlo Markov Chains (Dunkley et al. 2008). We use an exact maximum likelihood form for 
the l < 32 TT data (Eriksen et al. 2007b). We have investigated theoretically optimal methods for 
incorporating window function uncertainties into the likelihood, but in tests with simulated data we 
have found them to be biased. In the end, we adopt the form used in the 3 year analysis (Hinshaw 
et al. 2007), but we incorporate the smaller 5 year window function uncertainties (Hill et al. 2008) as 
inputs. We now routinely account for gravitational lensing when assessing parameters, and we have 
added an option to use low-/ TB and EB data for testing non-standard cosmological models. . 

• For testing nongaussianity, we employ an improved estimator for f^L (Creminelli et al. 2006; Yadav 
et al. 2007). The results of this analysis are described in Komatsu et al. (2008). 


3. OBSERVATIONS AND MAPS 

The 5 year WMAP data encompass the period from 00:00:00 UT, 10 August 2001 (day number 222) to 
00:00:00 UT, 9 August 2006 (day number 222). The observing efficiency during this time is roughly 99%; 
Table 2 lists the fraction of data that was lost or rejected as unusable. The Table also gives the fraction of 
data that is flagged due to potential contamination by thermal emission from Mars, Jupiter, Saturn, Uranus, 
and Neptune. These data are not used in map-making, but are useful for in-flight beam mapping (Hill et al. 
2008; Limon et al. 2008). 

After performing an end-to-end analysis of the instrument calibration, single-year sky maps are created 
from the time-ordered data using the procedure described by Jarosik et al. (2007). Figure 1 shows the 5 year 
temperature maps at each of the five WMAP observing frequencies: 23, 33, 41, 61, and 94 GHz. The number 
of independent observations per pixel, iVobs, is qualitatively the same as Figure 2 of Hinshaw et al. (2007) 
and is not reproduced here. The noise per pixel, p, is given by a(p) = &oN~^J 2 (p) , where <to is the noise 
per observation, given in Table 1. To a very good approximation, the noise per pixel in the 5 year maps is 
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a factor of y/E times lower than in the single-year maps. Figures 2 and 3 show the 5 year polarization maps 
in the form of the Stokes parameters Q and U, respectively. Maps of the relative polarization sensitivity, 
the Q and U analogs of A 0 b s , are shown in Figure 13 of Jarosik et al. (2007) and are not updated here. 
A description of the low-resolution pixel-pixel inverse covariance matrices used in the polarization analysis 
is also given in Jarosik et al. (2007), and is not repeated here. The polarization maps are dominated by 
foreground emission, primarily synchrotron emission from the Milky Way. Figure 4 shows the polarization 
maps in a form in which the color scale represents polarized intensity, P — \/Q 2 -\ -U 2 , and the line segments 
indicate polarization direction for pixels with a signal-to- noise ratio greater than 1. As with the temperature 
maps, the noise per pixel in the 5 year polarization maps is y/E times lower than in the single-year maps. 

Figure 5 shows the difference between the 5 year temperature maps and the corresponding 3 year maps. 
All maps have been smoothed to 2° resolution to minimize the noise difference between them (due to the 
additional years of data). The left column shows the difference without any further processing, save for the 
subtraction of a relative offset between the maps. Table 3 gives the value of the relative offset in each band. 
Recall that WMAP is insensitive to absolute temperature, so we adopt a convention that sets the zero level 
in each map based on a model of the foreground emission at the galactic poles. While we have not changed 
conventions, our 3 year estimate was erroneous due to the use of a preliminary CMB signal map at the time 
the estimate was made. This error did not affect any cosmological results, but it probably explains the offset 
differences noted by Eriksen et al. (2007a) in their recent analysis of the 3 year data. 

The dominant structure in the left column of Figure 5 consists of a residual dipole and galactic plane 
emission. This reflects the updated 5 year calibration which has produced changes in the gain of order 0.3% 
compared to the 3 year gain estimate (see §4 for a more detailed discussion of the calibration). Table 3 gives 
the dipole amplitude difference in each band, along with the much smaller quadrupole and octupole power 
difference. (For comparison, we estimate the CMB power at l = 2,3 to be /(/ -F l)Ci/2n = 211, 1041 /iK 2 , 
respectively.) The right column of Figure 5 shows the corresponding sky map differences after the 3 year 
map has been rescaled by a single factor (in each band) to account for the mean gain change between the 3 
and 5 year calibration determinations. The residual galactic plane structure in these maps is less than 0.2% 
of the nominal signal in Q band, and less than 0.1% in all the other bands. The large scale structure in the 
band- averaged temperature maps is quite robust. 


4. CALIBRATION IMPROVEMENTS 

With the 5 year processing we have refined our procedure for evaluating the instrument calibration, and 
have improved our estimates for the calibration uncertainty. The fundamental calibration source is still the 
dipole anisotropy induced by WMAP' s motion with respect to the CMB rest frame (Hinshaw et al. 2003; 
Jarosik et al. 2007), but several details of the calibration fitting have been modified. The new calibration 
solution is consistent with previous results in the overlapping time range. We estimate the uncertainty in 
the absolute calibration is now 0.2% per differencing assembly. 

The basic calibration procedure posits that a single channel of time-ordered data, d*, may be modeled 
as 

di — g% [A T v i -f- A T a i] bi, (1) 

where i is a time index, gi and b{ are the instrument gain and baseline, at time step z, A T v i is the differential 
dipole anisotropy induced by WMAP s motion, and AT^ is the differential sky anisotropy. We assume that 
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A T v i is known exactly and has the form 


A T vi = 


To 

c 


v» • [(1 +x im )nA,i - 


(1 — X\ m )liB,i], 


(2) 


where To = 2.725 K is the CMB temperature (Mather et al. 1999), c is the speed of light, v* is WMAP' s 
velocity with respect to the CMB rest frame at time step z, x- im is the loss imbalance parameter (Jarosik 
et al. 2007), and n^,i, and ns,i are the unit vectors of the A- and B-side lines of sight at time step i (in the 
same frame as the velocity vector). The velocity may be decomposed as 


Vi = VWMAP— SSB,i + VSSB-CMB, (3) 

where the first term is WMAP ' s velocity with respect to the solar system barycenter, and the second is the 
barycenter velocity with respect to the CMB. The former is well determined from ephemeris data, while the 
latter has been measured by COBE-DMR with an uncertainty of 0.7% (Kogut et al. 1996). Since the latter 
velocity is constant over WMAP ' s life span, any error in our assumed value of vssb-cmb will, in theory, 
be absorbed into a dipole contribution to the anisotropy map, T a . We test this hypothesis below. The 
differential sky signal has the form 

AT a j = (1 + %im)[I a (PA,i) T P a{PA,ii 7+i)] — (1 — %im)[Ia(PB,i) Ta(jPB,ii 7 B,i)]? (4) 

where pa,i is the pixel observed by the A-side at time step i (and similarly for B), I a (p) is the temperature 
anisotropy in pixel p (the intensity Stokes parameter, /), and P a (p , 7) is the polarization anisotropy in pixel 
p at polarization angle 7 (Hinshaw et al. 2003) which is related to the linear Stokes parameters Q and U by 


Pa(pn) = Q(p) cos 27 + 17 (p) sin 2y. (5) 

We further note that, in general, I a and P a depend on frequency owing to Galactic emission. 

A main goal of the data processing is to simultaneously fit for the calibration and sky signal. Unfor- 
tunately, since the data model is nonlinear and the number of parameters is large, the general problem is 
intractable. In practice, we proceed iteratively as follows. Initially we assume the gain and baseline are 
constant for a given time interval, typically between 1 and 24 hours, 


9i — Gk Tfc < t{ < T k . j_i (6) 

bi = B k T k <ti<T k + 1, (7) 


where t{ is the time of the zth individual observation, and r k is the start time of the kth calibration interval. 
Throughout the fit we fix the velocity-induced signal, equation (2), using vssb-cmb = [—26.29, —244.96, +275.93] 
km s -1 (in Galactic coordinates), and, for the first iteration, we assume no anisotropy signal, AT a = 0. Then, 
for each calibration interval k we perform a linear fit for G k and B k with fixed AT*, + AT a . As we pro- 
ceed through the intervals, we apply this calibration to the raw data and accumulate a new estimate of the 
anisotropy map as per equation 19 of Hinshaw et al. (2003). The procedure is repeated with each updated 
estimate of AT a . Once the calibration solution has converged, we fit the gain data, G k , to a model that 
is parameterized by the instrument detector voltage and the temperatures of the receiver’s warm and cold 
stages, equation 2 of Jarosik et al. (2007). This parametrization still provides a good fit, to the G k data, so 
we have not updated its form for the 5 year analysis. The updated best-fit parameters are given in the 5 
year Explanatory Supplement (Limon et al. 2008). Note that for each radiometer, the relative gain vs. time 
over 5 years is determined by just two parameters. 

For the 5 year processing we have focused on the veracity of the “raw” calibration, G k and B k . Specif- 
ically, we have improved and/or critically reexamined several aspects of the iterative fitting procedure: 
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• We have incorporated the effect of far sidelobe pickup directly into the iterative calibration procedure, 
rather than as a fixed correction (Jarosik et al. 2007). We do this by segregating the differential signal 
into a main beam contribution and a sidelobe contribution, 

AT* = AT ma i ni i + AT s id e ,i- (8) 

(Hill et al. 2008 discuss how this segregation is defined in the 5 year processing.) After each iteration 
of the calibration and sky map estimation, we (re)compute a database of AX^e on a grid of pointings 
using the new estimate of I a . We then interpolate the database to estimate AX s id e ,i for each time 
step i. Note that AT s id e includes contributions from both the velocity-induced signal and the intrinsic 
anisotropy. Ignoring sidelobe pickup can induce gain errors of up to 1.5% in K band, 0.4% in Ka band, 
and ~0.25% in Q-W bands. 

• In general, the different channels within a DA have different center frequencies (Jarosik et al. 2003); 
hence the different channels measure a slightly different anisotropy signal due to differences in the 
Galactic signal. We assess the importance of accounting for this in the calibration procedure. 

• A single DA channel is only sensitive to a single linear polarization state. ( WMAP measures polariza- 
tion by differencing orthogonal polarization channels.) Thus we cannot reliably solve for both P a and 
for I a at each channel’s center frequency. We assess the relative importance of accounting for one or 
the other on both the gain and baseline solutions. 

• We examine the sensitivity of the calibration solution to the choice of vssb-cmb and to assumptions 
of time-dependence in the gain. 


4.1. Calibration Tests 

We use a variety of end-to-end simulations to assess and control the systematic effects noted above. We 
summarize a number of the key tests in the remainder of this section. 

The first case we consider is a noiseless simulation in which we generate time-ordered data from an input 
anisotropy map which includes CMB and Galactic foreground signal (one map per channel, evaluated at the 
center frequency of each channel) and a known dipole amplitude. The input gain for each channel is fixed 
to be constant in time. We run the iterative calibration and sky map solver allowing for an independent sky 
map solution at each channel (but no polarization signal). When fitting for the calibration, we assume that 
vssb-cmb differs from the input value by 1% to see if the known, modulated velocity term, vwmap-ssb, 
properly “anchors” the absolute gain solution. The results are shown in the top panel of Figure 6 where it is 
shown that the absolute gain recovery is robust to errors in vssb-cmb- We recover the input gain to better 
than 0.1% in this instance. 

The second case we consider is again a noiseless simulation that now includes only dipole signal (with 
Earth- velocity modulation), but here we vary the input gain using the flight-derived gain model (Jarosik 
et al. 2007). The iterative solver was run on the K band data for 1400 iterations, again starting with an 
initial guess that was in error by 1%. The results are shown in the bottom panel of Figure 6, which indicate 
systematic convergence errors of >0.3% in the fitted amplitude of the recovered gain model. Since the 
input sky signal in this case does not have any Galactic foreground or polarization components, we cannot 
ascribe the recovery errors to the improper handling of those effects in the iterative solver. We have also 
run numerous other simulations that included various combinations of instrument noise, CMB anisotropy, 
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Galactic foreground signal (with or without individual center frequencies per channel), polarization signal, 
and input gain variations. The combination of runs are too numerous to report on in detail, and the results 
are not especially enlightening. The most pertinent trend we can identify is that when the input value of 
v ssb-cmb is assumed in the iterative solver, the recovered gain is in good agreement with the input, but 
when the initial guess is in error by 1%, the recovered gain will have comparable errors. We believe the 
lack of convergence is due to a weak degeneracy between gain variations and the sky map solution. Such a 
degeneracy is difficult to diagnose in the context of this iterative solver, especially given the computational 
demands of the system, so we are assessing the system more directly with a low-resolution parameterization 
of the gain and sky signal, as outlined in Appendix A. 

Since the latter effort is still underway, we have adopted a more pragmatic approach to evaluating the 
absolute gain and its uncertainty for the 5 year data release. We proceed as follows: after 50 iterations of the 
calibration and sky map solver, the dominant errors in the gain and sky map solution are 1) a dipole in the 
sky map, and 2) a characteristic wave form that reflects a relative error between vssb-cmb and vwmap-ssb- 
At this point we can calibrate the amplitude of the gain error wave form to the magnitude of the velocity 
error in vssb-cmb- We can then fit the gain solution to a linear combination of the gain model of Jarosik 
et al. (2007) and the velocity error wave form. See Appendix B for details on this fitting procedure. In 
practice this fit is performed simultaneously on both channels of a radiometer since those channels share 
one gain model parameter. We have tested this procedure on a complete flight-like simulations that includes 
every important effect known, including input gain variations. The results of the gain recovery are shown in 
Figure 7, and based on this we conservatively assign an absolute calibration uncertainty of 0.2% per channel 
for the 5 year WMAP archive. 


4.2. Summary 

The series of steps taken to arrive at the final 5 year calibration are as follows: 

• Run the iterative calibration and sky map solver over the full 5 year data set for 50 iterations, using 
24 hour calibration intervals. This run starts with I a = P a = 0 and updates I a for each individual 
channel of data. P a is assumed to be 0 throughout this run. We keep the gain solution, Gk, from this 
run and discard the baseline solution. 

• Run the iterative calibration and sky map solver over the full 5 year data set for 50 iterations, using 
1 hour calibration intervals. This run starts with I a = P a = 0 and updates both using the intensity 
and polarization data in the two radiometers per DA, as per Appendix D of Hinshaw et al. (2003). 
We keep the baseline solution, Bk , from this run and discard the gain solution. Both of these runs 
incorporate the sidelobe correction as noted above. 

• Fit the gain solution, Gk simultaneously for the gain model and for an error in the velocity, Avssb-cmb, 
as described in Appendix B. This fit is performed on two channels per radiometer with the gain model 
parameter To common to both channels. 

• We average the best-fit velocity error over all channels within a frequency band under the assumption 
that the dipole is the same in each of these channels. We then fix the velocity error to a single value 
per frequency band and re-fit the gain model parameters for each pair of radiometer channels. 
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Based on end-to-end simulations with flight-like noise, we estimate the absolute gain error per radiometer 
to be 0.2%. We believe the limiting factor in this estimate is a weak degeneracy between thermal variations 
in the instrument gain, which are annually modulated, and annual variations induced by errors in vssb-cmb- 
Since there is a small monotonic increase in the spacecraft temperature, additional years of data should allow 
improvements in our ability to separate these effects. 

Once we have finalized the gain model, we form a calibrated time-ordered data archive using the gain 
model and the 1 hour baseline estimates to calibrate the data. This archive also has a final estimate of the 
far sidelobe pickup subtracted from each time-ordered data point. However, we opt not to subtract a dipole 
estimate from the archive at this stage in the processing. 


5. BEAM IMPROVEMENTS 

In addition to reassessing the calibration, the other major effort undertaken to improve the 5 year data 
processing was to extend the physical optics model of the WMAP telescope based on flight measurements of 
Jupiter. This work is described in detail in Hill et al. (2008) so we only summarize the key results with an 
emphasis on their scientific implications. The basic aim of the work is to use the flight beam maps from all 
10 DA’s to determine the in-flight distortion of the mirrors. This program was begun for the A-side mirror 
during the 3 year analysis; for the 5 year analysis we have quadrupled the number of distortion modes we 
fit (probing distortion scales that are half the previous size), and we have developed a completely new and 
independent model of the B-side distortions, rather than assuming that they mirror the A-side distortions. 
We have also placed limits on smaller scale distortions by comparing the predicted beam response at large 
angles to sidelobe data collected during WMAP 1 s early observations of the Moon. 

Given the best-fit mirror model, we compute the model beam response for each DA and use it in 
conjunction with the flight data to constrain the faint tails of the beams, beyond ~ 1° from the beam peak. 
These tails are difficult to constrain with flight data alone because the Jupiter signal to noise ratio is low, but, 
due to their large areal extent they contain a non-negligible fraction (up to 1%) of the total beam solid angle. 
An accurate determination of the beam tail is required to properly measure the ratio of sub- degree- scale 
power to larger-scale power in the diffuse CMB emission (and to accurately assign point source flux). 

Figure 14 in Hill et al. (2008) compares the beam radial profiles used in the 3 year and 5 year analyses, 
while Figure 13 compares the /-space transfer functions derived from the Legendre transform of the radial 
profile. The important changes to note are the following. 

1. In both analyses we split the beam response into main beam and far sidelobe contributions. In the 5 
year analysis we have enlarged the radius at which this transition is made (Hill et al. 2008). In both 
cases, we correct the time-ordered data for far sidelobe pickup prior to making sky maps, while the 
main beam contribution is only accounted for in the analysis of sky maps, e.g., in power spectrum 
deconvolution. As a result, the sky maps have a slightly different effective resolution which is most 
apparent in K band, as in Figure 5. However, in each analysis, the derived transfer functions are 
appropriate for the corresponding sky maps. 

2. In the 3 year analysis, the main beam profile was described by a Hermite polynomial expansion fit to 
the observations of Jupiter in the time-ordered data. This approach was numerically problematic in 
the 5 year analysis due to the larger transition radius; as a result, we now simply co-add the time- 
ordered data into radial bins to obtain the profiles. In both cases, the underlying time-ordered data 
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is a hybrid archive consisting of flight data for points where the beam model predicts a value above 
a given contour, and model values for points below the contour (Hill et al. 2008). With the improved 
beam models and a new error analysis, we have adjusted these hybrid contours down slightly, with the 
result that we use proportionately more flight data (per year) in the new analysis. The radius at which 
the 5 year profile becomes model dominated (>50% of the points in a bin) is indicated by dotted lines 
in Figure 14 of Hill et al. (2008). 

3. The right column of Figure 14 in Hill et al. (2008) shows the fractional change in solid angle due to 
the updated profiles. The main point to note is the ~1% increase in the V2 and W band channels, 
primarily arising in the bin from 1 to 2 degrees off the beam peak. As can be seen in Figure 3 of Hill 
et al. (2008), this is the angular range in which the new beam models produced the most change, owing 
to the incorporation of smaller distortion modes in the mirror model. The 3 year analysis made use 
of the model in this angular range which, in hindsight, was suppressing up to ~1% of the solid angle 
in the V and W band beams. (The longer wavelength channels are less sensitive to distortions in this 
range, so the change in solid angle is smaller for K-Q bands.) In the 5 year analysis, we use relatively 
more flight data in this regime, so we are less sensitive to any remaining model uncertainties. Hill 
et al. (2008) place limits on residual model errors and propagate those errors into the overall beam 
uncertainty. 

4. Figure 13 in Hill et al. (2008) compares the beam transfer functions, derived by transforming the 3 
year and 5 year radial profiles. (To factor out the effect of changing the transition radius, the 3 year 
profiles were extended to the 5 year radius using the far sidelobe data, for this comparison.) Since the 
transfer functions are normalized to 1 at / = 1, the change is restricted to high /. In V and W bands, 
bi has decreased by ~0.5 - 1% due largely to the additional solid angle picked up in the 1-2 degree 
range. This amounts to a ~1 a change in the functions, as indicated by the red curves in the Figure. 

The calibrated angular power spectrum is proportional to l/<7 2 6 2 , where g is the mean gain and bi is the 
beam transfer function, thus the net effect of the change in gain and beam determinations is to increase the 
power spectrum by ~0.5% at l < 100, and by ~2.5% at high /. Nolta et al. (2008) give a detailed evaluation 
of the power spectrum while Dunkley et al. (2008) and Komatsu et al. (2008) discuss the implications for 
cosmology. 


6. LOW -l POLARIZATION TESTS 

The 3 year data release included the first measurement of microwave polarization over the full sky, in the 
form of Stokes Q and U maps in each of 5 bands. The analysis of WMAP polarization data is complicated by 
the fact that the instrument was not designed to be a true polarimeter, thus a number of systematic effects 
had to be understood prior to assigning reliable error estimates to the data. Page et al. (2007) presented 
the 3 year polarization data in great detail. In this section we extend that analysis by considering some 
additional tests that were not covered in the 3 year analysis. 


6.1. Year-to-Year Consistency Tests 

With 5 years of data it is now possible to subject the data to more stringent consistency tests than was 
previously possible. In general, the number of independent cross-power spectra we can form within a band 
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with N d differencing assemblies is N d (N d - l)/2 x N y + Nj x N y (N y - l)/2. With 5 years of data, this gives 
10 independent estimates each in K and Ka band, 45 each in Q and V band, and 190 in W band. For cross 
power spectra of distinct band pairs, with N d i and N d 2 DA’s in each band, the number is N d iN d 2 x N y . This 
gives 50 each in KaQ and KaV, 100 each in KaW and QV, and 200 each in QW and VW. (For comparison, 
the corresponding numbers are 3, 15, & 66, and 18, 36, & 72 with 3 years of data.) 

We have evaluated these individual spectra from the 5 year data and have assigned noise uncertainties 
to each estimate using the Fisher formalism described in Page et al. (2007). We subject the ensemble to an 
internal consistency test by computing the reduced x 2 of the data at each multipole / within each band or 
band pair, under the hypothesis that the data at each multipole and band measures the same number from 
DA to DA and year to year. The results of this test are given in Table 4 for the foreground-cleaned EE, EB, 
and BB spectra from l = 2 — 10 for all band pairs from KaKa to WW. There are several points to note in 
these results. 

1. For l > 6, the most significant deviation from 1 in reduced x 2 > in an y spectrum or band, is 1.594 in 
the l = 7 BB spectrum for KaQ. With 50 degrees of freedom, this is a 3 a deviation, but given that 
we have 150 /> 6 samples in the table, we expect of order 1 such value. Thus we conclude that the 
Fisher- based errors provide a good description of the DA-to-DA and year-to-year scatter in the l > 6 
polarization data. If anything, there is a slight tendency to overestimate the uncertainties at higher /. 

2. For l < 5, we find 37 out of 120 points where the reduced x 2 deviates from 1 at more than 4 a 
significance, indicating excessive internal scatter in the data relative to the Fisher errors. However, 
all but 5 of these occur in cross-power spectra in which one or both of the bands contain W band 
data. If we exclude combinations with W band, the remaining 72 points have a mode in the reduced 
X 2 distribution of 1 with a slight positive skewness due to the 5 points noted above, which all contain 
Q band data. This may be a sign of slight foreground residuals contributing additional noise to the 
Q band data, though we do not see similar evidence in the Ka band spectra which would be more 
foreground contaminated prior to cleaning. For Ka-V bands, we believe that the Fisher errors provide 
an adequate description of the scatter in this l < 5 polarization data, but we subject polarization 
sensitive cosmological parameter estimates, e.g., the optical depth, to additional scrutiny in §6.3. 

3. Of special note is l = 3 BB which, as noted in Page et al. (2007), is the power spectrum mode that 
is least modulated in the WMAP time-ordered data. This mode is therefore quite sensitive to how 
the instrument baseline is estimated and removed and, in turn, to how the 1/f noise is modeled. In 
the accounting above, the l = 3 BB data have the highest internal scatter of any low-/ polarization 
mode. In particular, every combination that includes W band data is significantly discrepant; and the 
two most discrepant non-W band points are also estimates of l = 3 BB. We comment on the W band 
data further below, but note here that the final co-added BB spectrum (based on Ka, Q, and V band 
data) does not lead to a significant detection of tensor modes. However, we caution that any surprising 
scientific conclusions which rely heavily on the WMAP l = 3 BB data should be treated with caution. 

Based on the analysis presented above, we find the W band polarization data is still too unstable at 
low-/ to be reliably used for cosmological studies. We cite more specific phenomenology and consider some 
possible explanations in the remainder of this section. 

The 5 year co-added W band EE spectrum is shown in Figures 8, in the form of likelihood profiles 
from 1 — 2 — 1. At each multipole we show two curves: an estimate based on evaluating the likelihood 
multipole by multipole, and an estimate based on the pseudo-Q method (Page et al. 2007). The best-fit 
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model EE spectrum, based on the combined Ka, Q, and V band data is indicated by the dashed lines in each 
panel. Both spectrum estimates show excess power relative to the model spectrum, with the most puzzling 
multipole being l = 7 which, as shown in Table 4, has an internal reduced y 2 of 1.015, for 190 degrees of 
freedom. This data has the hallmark of a sky signal, but that hypothesis is implausible for a variety of 
reasons (Page et al. 2007). It is more likely due to a systematic effect that is common to a majority of the 
W band channels over a majority of the 5 years of data. We explore and rule out one previously neglected 
effect in §6.2. It is worth recalling that l = 7 EE, like l = 3 BB, is a mode that is relatively poorly measured 
by WMAP, as discussed in Page et al. (2007); see especially Figure 16 and its related discussion. 

The W band BB data also exhibit unusual behavior at l = 2, 3. In this case, these two multipoles have 
internal reduced y 2 greater than 6, and the co-added l = 2 point is nearly 10 a from zero. However, with 
190 points in each 5 year co-added estimate it is now possible to look for trends within the data that were 
relatively obscure with only 3 years of data. In particular, we note that in the l = 2 estimate, there are 28 
points that are individually more than 5 a from zero and that all of them contain W 1 data in one or both of 
the DA pairs in the cross power spectrum. Similarly for l = 3, there are 14 points greater than 5 a and all 
of those points contain W4 data in one or both of the DA pairs. We have yet to pinpoint the significance of 
this result, but we plan to study the noise properties of these DA’s beyond what has been reported to date, 
and to sharpen the phenomenology with additional years of data. 


6.2. Emissivity Tests 

In this section we consider time dependent emission from the WMAP optics as a candidate for explaining 
the excess W band “signal” seen in the EE spectrum, mostly at l = 7. In the end, the effect proved not to 
be significant, but it provides a useful illustration of a common-mode effect that we believe is still present 
in the W band polarization data. 

From a number of lines of reasoning, we know that the microwave emissivity of the mirrors is a few 
percent in W band, and that it scales with frequency roughly like u 15 across the WMAP frequency range, 
as expected for a classical metal (Born & Wolf 1980). Hence this mechanism has the potential to explain 
a common-mode effect that is primarily seen in W band. Further, Figure 1 in Jarosik et al. (2007) shows 
that the physical temperature of the primary mirrors are modulated at the spin period by ~200 K, with a 
dependence on solar azimuth angle that is highly repeatable from year to year. We believe this modulation is 
driven by solar radiation diffracting around the WMAP sun shield reaching the tops of the primary mirrors, 
which are only a few degrees within the geometric shadow of the sun shield. In contrast, the secondary 
mirrors and feed horns are in deep shadow and show no measurable variation at the spin period, so that any 
emission they produce only contributes to an overall radiometer offset, and will not be further considered 
here. 


As a rough estimate, the spin modulated emission from the primary mirrors could produce as much 
as ~ 0.02 x 200 = 4 fiK of radiometric response in W band, but the actual signal depends on the relative 
phase of the A and B-side mirror variations and the polarization state of the emission. In more detail, the 
differential signal, d(t), measured by a radiometer with lossy elements is 

d(t) = (1 - e A ) T A (t) - (1 - e B ) T B (t) + e p A T p (t ) - e p B T p (t) (9) 

where ea = e p A + e s A + e f A is the combined loss in the A-side optics: (p)rimary plus (s)econdary mirrors, plus 
the (f)eed horn, and likewise for the B-side. Ta,b is the sky temperature in the direction of the A or B-side 
line-of-sight; and T A B is the physical temperature of the A or B-side primary mirror. 
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The first two terms are the sky signal attenuated by the overall loss in the A and B side optics, 
respectively. The effects of loss imbalance, which arise when ea ^ £b, have been studied extensively (Jarosik 
et al. 2003, 2007). We account for loss imbalance in the data processing and we marginalize over residual 
uncertainties in the imbalance coefficients when we form the pixel-pixel inverse covariance matrices (Jarosik 
et al. 2007). Updated estimates of the loss imbalance coefficients based on fits to the 5 year data are reported 
in Table 5. 


In the remainder of this section we focus on the last two emissive terms in Equation 9. Recall that a 
WMAP differencing assembly consists of two radiometers, 1 and 2, that are sensitive to orthogonal linear po- 
larization modes. The temperature and polarization signals are extracted by forming the sum and difference 
of the two radiometer outputs; thus, the emission terms we need to evaluate are 


dw ±d p 2 (t) = n - n 


(10) 


where e p Al is the A-side primary mirror emissivity measured by radiometer 1, and so forth. The factor of 
1 — e in the denominator applies a small correction for the mean loss, e = (ea H- eb)/2, and arises from the 
process of calibrating the data to a known sky brightness temperature (§4). Note that we only pick up a 
polarized response if e\ ^ C2- 

We have simulated this signal in the time-ordered data using the measured primary mirror temperatures 
as template inputs. The emissivity coefficients were initially chosen to be consistent with the loss imbalance 
constraints. However, in order to produce a measurable polarization signal, we had to boost the emissivity 
differences to the point where they became unphysical, that is |ei —62! > |ei + €2 1- Nonetheless, it was 
instructive to analyze this simulation by binning the resulting data (which also includes sky signal and 
noise) as a function of solar azimuth. The results are shown in the top panel of Figure 9 which shows 3 years 
of co-added W band polarization data, the d\ — c/2 channel; the input emissive signal is shown in red for 
comparison. We are clearly able to detect such a signal with this manner of binning. We also computed the 
low-/ polarization spectra and found that, despite the large spin modulated input signal, the signal induced 
in the power spectrum was less than 2 fiK 2 in /(/ + 1) Cf E /2 k, which is insufficient to explain the / = 7 
feature in the W band EE spectrum. 

In parallel with the simulation analysis, we have binned the flight radiometer data by solar azimuth 
angle to search for spin modulated features in the polarization data. The results for W band are shown in the 
bottom panel of Figure 9 for the 5 year data. While the \ 2 P er degree of freedom relative to zero is slightly 
high, there is no compelling evidence for a coherent spin modulated signal at the ~2 11 K level. In contrast, 
the simulation yielded spin modulated signals of 5-10 /zK and still failed to produce a significant effect in the 
EE spectrum. Hence we conclude that thermal emission from the WMAP optics cannot explain the excess 
W band EE signal. In any event, we continue to monitor the spin modulated data for the emergence of a 
coherent signal. 


6.3. Ka Band Tests 

The analysis presented in §6.1 shows that the Ka band polarization data is comparable to the Q and V 
band data in its internal consistency. That analysis was performed on data that had been foreground cleaned 
using the template method discussed in Page et al. (2007) and updated in Gold et al. (2008). In order to 
assess whether or not this cleaned Ka band data is suitable for use in cosmological parameter estimation we 
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subject it to two further tests: 1) a null test in which Ka band data is compared to the combined Q and V 
band data, and 2) a parameter estimation based solely on Ka band data. 

For the null test, we form polarization maps by taking differences, |Ska — where S = Q,U are 

the polarization Stokes parameters, Ska are the maps formed from the Ka band data, and Sqv are the maps 
formed from the optimal combination of the Q and V band data. We evaluate the EE power spectrum from 
these null maps by evaluating the likelihood mode by mode while holding the other multipoles fixed at zero. 
The results are shown in Figure 10, along with the best-fit model spectrum based on the final 5 year ACDM 
analysis. The spectrum is clearly consistent with zero, but to get a better sense of the power of this test, 
we have also used these null maps to estimate the optical depth parameter, r. The result of that analysis is 
shown as the dashed curve in Figure 11, where we find that the null likelihood peaks at r = 0 and excludes 
the most-likely cosmological value with ~95% confidence. 

As a separate test, we evaluate the r likelihood using only the template-cleaned Ka band signal maps. 
The result of that test is shown as the blue curve in Figure 11. While the uncertainty in the Ka band estimate 
is considerably larger than the combined QV estimate (shown in red), the estimates are highly consistent. 
The result of combining Ka, Q, and V band data is shown in the black curve. 

Dunkley et al. (2008) present a complementary method of foreground cleaning that makes use of Ka 
band data, in conjunction with K, Q, and V band data. Using a full 6 parameter likelihood evalutaion, 
they compare the optical depth inferred from the two cleaning methods while using the full combined data 
sets in both cases: see Figure 9 of Dunkley et al. (2008) for details. Based on these tests, we conclude that 
the Ka band data is sufficiently free of systematic errors and residual foreground signals that it is suitable 
for cosmological studies. The use of this band significantly enhances the overall polarization sensitivity of 
WMAP. 


7. SUMMARY OF 5- YEAR SCIENCE RESULTS 

Detailed presentations of the scientific results from the 5 year data are given by Gold et al. (2008), 
Wright et al. (2008), Nolta et al. (2008), Dunkley et al. (2008), and Komatsu et al. (2008). Starting with 
the 5 year temperature and polarization maps, with their improved calibration, Gold et al. (2008) give a 
new Markov Chain Monte Carlo-based analysis of foreground emission in the data. Their results are broadly 
consistent with previous analyses by the WMAP team and others (Eriksen et al. 2007b), while providing 
some new results on the microwave spectra of bright sources in the Galactic plane that aren’t well fit by 
simple power-law foreground models. Figure 12 shows the 5 year CMB map based on the internal linear 
combination (ILC) method of foreground removal. 

Wright et al. (2008) give a comprehensive analysis of the extragalactic sources in the 5 year data, 
including a new analysis of variability made possible by the multi-year coverage. The 5 year WMAP source 
catalog now contains 390 objects and is reasonably complete to a flux of 1 Jy away from the Galactic plane. 
The new analysis of the WMAP beam response (Hill et al. 2008) has led to more precise estimates of the 
point source flux scale for all 5 WMAP frequency bands. This information is incorporated in the new source 
catalog (Wright et al. 2008), and is also used to provide new brightness estimates of Mars, Jupiter, and 
Saturn (Hill et al. 2008). We find significant (and expected) variability in Mars and Saturn over the course 
of 5 years and use that information to provide a preliminary recalibration of a Mars brightness model (Wright 
2007), and to fit a simple model of Saturn’s brightness as a function of ring inclination. 
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The temperature and polarization power spectra are presented in Nolta et al. (2008). The spectra are 
all consistent with the 3 year results with improvements in sensitivity commensurate with the additional 
integration time. Further improvements in our understanding of the absolute calibration and beam response 
have allowed us to place tighter uncertainties on the power spectra, over and above the reductions from 
additional data. These changes are all reflected in the new version of the WMAP likelihood code. The most 
notable improvements arise in the third acoustic peak of the TT spectrum, and in all of the polarization 
spectra; for example, we now see unambiguous evidence for a 2nd dip in the high-/ TE spectrum, which 
further constrains deviations from the standard ACDM model. The 5 year TT and TE spectra are shown 
in Figure 13. We have also generated new maximum likelihood estimates of the low-/ polarization spectra: 
TE, TB, EE, EB, and BB to complement our earlier estimates based on pseudo -Ci methods (Nolta et al. 
2008). The TB, EB, and BB spectra remain consistent with zero. 

The cosmological implications of the 5 year WMAP data are discussed in detail in Dunkley et al. (2008) 
and Komatsu et al. (2008). The now-standard cosmological model: a flat universe dominated by vacuum 
energy and dark matter, seeded by nearly scale- invariant, adiabatic, Gaussian random-phase fluctuations, 
continues to fit the 5 year data. WMAP has now determined the key parameters of this model to high 
precision; a summary of the 5 year parameter results is given in Table 6. The most notable improvements 
are the measurements of the dark matter density, f 2 c /i 2 , and the amplitude of matter fluctuations today, 
erg. The former is determined with 6% uncertainty using WMAP data only (Dunkley et al. 2008), and with 
3% uncertainty when WMAP data is combined with BAO and SNe constraints (Komatsu et al. 2008). The 
latter is measured to 5% with WMAP data, and to 3% when combined with other data. The redshift of 
reionization is z r eion = 11.0 ± 1.4, if the universe were reionized instantaneously. The 2 cr lower limit is 
2reion>8-2, and instantaneous reionization at z re \ on = 6 is rejected at 3.5 a. The WMAP data continues to 
favor models with a tilted primordial spectrum, n s — 0.963^0 oi5- Dunkley et al. (2008) discuss how the 
ACDM model continues to fit a host of other astronomical data as well. 

Moving beyond the standard ACDM model, when WMAP data is combined with BAO and SNe obser- 
vations (Komatsu et al. 2008), we find no evidence for running in the spectral index of scalar fluctuations, 
dn s /d\nk = — 0.032^o 020 (^8% CL). The new limit on the tensor-to-scalar ratio is r < 0.20 (95% CL), and 
we obtain tight, simultaneous limits on the (constant) dark energy equation of state and the spatial curvature 
of the universe: —0.11 < 1 + w < 0.14 (95% CL) and —0.0175 < D*, < 0.0085 (95% CL). The angular power 
spectrum now exhibits the signature of the cosmic neutrino background: the number of relativistic degrees 
of freedom, expressed in units of the effective number of neutrino species, is found to be AT e ff = 4.4 ± 1.5 (68% 
CL), consistent with the standard value of 3.04. Models with iV e ff = 0 are disfavored at >99.5% confidence. 
A summary of the key cosmological parameter values is given in Table 6, where we provide estimates using 
WMAP data alone and WMAP data combined with BAO and SNe observations. A complete tabulation of 
all parameter values for each model and dataset combination we studied is available on LAMBDA. 

The new data also place more stringent limits on deviations from Gaussianity, parity violations, and 
the amplitude of isocurvature fluctuations (Komatsu et al. 2008). For example, new limits on physically 
motivated primordial non-Gaussianity parameters are —9 < /Jy£ al < 111 (95% CL) and —151 < < 253 

(95% CL) for the local and equilateral models, respectively. 
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8. CONCLUSIONS 

We have presented an overview of the 5 year WMAP data and have highlighted the improvements we 
have made to the data processing and analysis since the 3 year results were presented. The most substantive 
improvements to the processing include a new method for establishing the absolute gain calibration (with 
reduced uncertainty), and a more complete analysis of the WMAP beam response made possible by additional 
data and a higher fidelity physical optics model. Numerous other processing changes are outlined in §2. 

The 5 year sky maps are consistent with the 3 year maps and have noise levels that are y/h times less 
than the single year maps. The new maps are compared to the 3 year maps in §3. The main changes to 
the angular power spectrum are as follows: at low multipoles (/ < 100) the spectrum is ~0.5% higher than 
the 3 year spectrum (in power units) due to the new absolute gain determination. At higher multipoles it is 
increased by ~2.5%, due to the new beam response profiles, as explained in §5 and in Hill et al. (2008). These 
changes are consistent with the 3 year uncertainties when one accounts for both the 0.5% gain uncertainty 
(in temperature units) and the 3 year beam uncertainties, which were incorporated into the likelihood code. 

We have applied a number of new tests to the polarization data to check internal consistency and to 
look for new systematic effects in the W band data (§6). As a result of these tests, and of new analyses 
of polarized foreground emission (Dunkley et al. 2008), we have concluded that Ka band data can be used 
along with Q and V band data for cosmological analyses. However, we still find a number of features in the 
W band polarization data that preclude its use, except in the Galactic plane where the signal to noise is 
relatively high. We continue to investigate the causes of this and have identified new clues to follow up on 
in future studies (§6.1). 

Scientific results gathered from the suite of 5 year papers are summarized in §7. The highlights include 
smaller uncertainties in the optical depth, r, due to a combination of additional years of data and to the 
inclusion of Ka band polarization data: instantaneous reionization at z Teion = 6 is now rejected at 3.5 cr. 
New evidence favoring a non-zero neutrino abundance at the epoch of last scattering, made possible by 
improved measurements of the third acoustic peak; and new limits on the nongaussian parameter Jnl, 
based on additional data and the application of a new, more optimal bispectrum estimator. The 5 year 
data continue to favor a tilted primordial fluctuation spectrum, in the range n s ~ 0.96, but a purely scale 
invariant spectrum cannot be ruled out at >3 a confidence. 

The WMAP observatory continues to operate at L2 as designed, and the addition of two years of 
flight data has allowed us to make significant advances in characterizing the instrument. Additional data 
beyond 5 years will give us a better understanding of the instrument, especially with regards to the W band 
polarization data since the number of jackknife combinations scales like the square of the number of years 
of operation. If W band data can be incorporated into the EE power spectrum estimate, it would become 
possible to constrain a second reionization parameter and thereby further probe this important epoch in 
cosmology. The WMAP data continues to uphold the standard ACDM model but more data may reveal 
new surprises. 


9. DATA PRODUCTS 

All of the WMAP data is released to the research community for further analysis through the Legacy 
Archive for Microwave Background Data Analysis (LAMBDA) at http://lambda.gsfc.nasa.gov. The products 
include the complete 5 year time-ordered data archive (both raw and calibrated); the calibrated sky maps 
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in a variety of processing stages (single year by DA, multi-year by band, high resolution and low resolution, 
smoothed, foreground-subtracted, and so forth); the angular power spectra and cosmological model likelihood 
code; a full table of model parameter values for a variety of model and data sets (including the best-fit 
model spectra and Markov chains); and a host of ancillary data to support further analysis. The WMAP 
Explanatory Supplement provides detailed information about the WMAP in-flight operations and data 
products (Limon et al. 2008). 


The WMAP mission is made possible by the support of the Science Mission Directorate Office at NASA 
Headquarters. This research was additionally supported by NASA grants NNG05GE76G, NNX07AL75G 
SOI, LTSA03-000-0090, ATPNNG04GK55G, and ADP03-0000-092. EK acknowledges support from an Al- 
fred P. Sloan Research Fellowship. This research has made use of NASA’s Astrophysics Data System Bibli- 
ographic Services. We acknowledge use of the HEALPix, CAMB, CMBFAST, and CosmoMC packages. 
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A. FISHER MATRIX ANALYSIS OF CALIBRATION AND SKY MAP FITS 

A.l. Least Squares Calibration and Sky Model Fitting 

Let i be a time index in the time ordered data. Let gi be parameters for the gain, ai m be parameters 
for the temperature anisotropy and b k be parameters for the baseline offset. 

The model of the time-ordered data (TOD) is 


wii — Qi [AT^j T AX^j] T 


(Al) 


where i is a time index, A T v { is the differential dipole signal at time step i, including the CMB dipole, and 
A T ai is the differential anisotropy signal at time step i. The parameters of the model are the hourly gain 
and baseline values, and the sky map pixel temperatures (which goes into forming AX a . We fit for them by 
minimizing 

( Ci - rrij ) 2 




=x 


(A2) 


where c* is the raw data, in counts, and cq is the rms of the ith observation, in counts. The Fisher matrix 
requires taking the second derivative of y 2 with respect to all parameters being fit. In order to reduce the 
dimensionality of the problem to something manageable, we expand the calibration and sky signal in terms 
of a small number of parameters. We can write 


j 

(A3) 

k 

(A4) 

^ ^ O’lm [^Zm(h'Ai) Pzm(f^Bi)] i 

Im 

(A5) 


This preprint was prepared with the A AS IAT]eX macros v5.2. 
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where G and B are function of time (defined below), aim are the harmonic coefficients of the map, and fiAi 
is the unit vector of the A-side feed at time step z, and likewise for B. 


A reasonable set of basis functions for the gain and baseline allow for an annual modulation and a small 
number of higher harmonics. Note that this does not include power at the spin or precession period, which 
might be an important extension to consider. For now we consider the trial set 


and 


Gji = 


1 

cos jOi 

s kl(j — Jmax)$i 


j = 0 
j = 1, ■ 


? Jmax 


J — Jmax 


+ 1..-..2 j 


max 


Bki = 


1 

cos kOi 

sin (/c ^max)^i 


k = 0 

k = 1 , . . . , k mSLX 

k — ^max + 1, • . . , *2k n 


(A6) 


(A7) 


where 0 = tan 1 (h y /h x ). Here n is the unit vector from WMAP to the Sun, and the components are 
evaluated in ecliptic coordinates. 


A. 2. Evaluation of the Fisher Matrix 


We wish to evaluate the 2nd derivative 

1 d 2 x 2 

2 dpidpj 

where pi and pj are the parameters we are trying to fit. The needed first derivatives are 

1 dx 2 _ ( C i ~ m i)Gj'i [A Tyi + A Tai] 

2 d g y ^ a? 


Then 


1 <9x 2 

2 db k ' 



■ rrii)B k 


1 <?X 2 = y- (cj - mj)gj - Yi> m >(h B i)] 

2 dai'm' of 

l 1 

1 d 2 x 2 sr- Gyi [A T*i + A Tgj\ G i"i [ + Ar ai ] 

2 dg J 'dgJ" -y erf 


1 ^ 2 x 2 

2 dgi' dai'm' 


V - "- 9i [Yl* m f (flAi) Yl' m' (Agj)] Gj'i [ AX^i -f- AT* a j] 


+ C>y^(ci - mi) 

i 


1 c> 2 x 2 y^ ByjGfj [ATyj + Ar o j] 

2 dgi' db k ' of 

1 d 2 X 2 _ y- ffi [Vj'm'(nAi) - Aj'm'(wgi)] ffi [fy'm"(nAt) ~ 

2 dai’m' dai" m " . erf 

l 1 


1 d X y' 9iB k ' i m / (flAi) Yl'm'iaBi)] 

2 dai'm’ db k ' <x? 

l 1 


(A8) 

(A9) 

(A10) 

(All) 

(A12) 

(A13) 

(A14) 

(A15) 

(A16) 
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■1 d 2 X 2 _ Bk"jBk’i 
2 db k> db k " 2s (j ? 

7. * 


(A17) 


From this we can form the inverse covariance matrix 


/ i a 2 x 2 

/ 2 dgj'dgi" 


1 

2 dgj' ddin^n 


1 \ 

2dgi'db k " \ 


C- 1 = 


1 ay , 1 aV 1 a 2 x 2 

2 da lfrn /dgi" 2 da t i m / daj// m // 2 da l / rn / db k>> 


V 


1 aV 

2 db k 'dgj" 


1 flV 

2 db k> da t n m n 


1 / 

2 db k> db k " / 


(A18) 


where the gain and baseline blocks are (2j max + 1) x (2j max + 1), and the sky map block is (/m ax + l) 2 x 

(^max + l) 2 - 


If we decompose C 1 using SVD the parameter covariance matrix can be inverted to have the form 

c = Y.~X i] ® V{i) {A19) 

i 

where the Wi are the singular values, and the are the columns of the orthogonal matrix V. In this form, 
the uncertainty in the linear combination of parameters defined by is l/wi. 


B. CALIBRATION MODEL FITTING WITH GAIN ERROR TEMPLATES 


B.l. Gain Error From Calibration Dipole Error 

Consider a simple model where the input sky consists of only a pure fixed (CMB) dipole, described by 
the vector d c , and a dipole modulated by the motion of WMAP with respect to the Sun, described by the 
time-dependent vector d i; (t). The raw data produced by an experiment observing this signal is 

c(U ) = g(ti)[At c (ti) + A t v (ti)\ (Bl) 


where c(ti) is the TOD signal in counts, g(ti) is the true gain of the instrument and is the differential 

signal produced by each dipole component (m = c, v ) at time ti given the instrument pointing at that time. 
Note that we have suppressed the explicit baseline and noise terms here for simplicity. 


Now suppose we calibrate the instrument using an erroneous CMB dipole, d' c — rd c = (1 + Ar)d c , 
where r is a number of order one (and Ar « 1 so we can ignore terms of order Ar 2 ). The fit gain, #/(£), 
will then roughly have the form 


9 fit) = 


eft) 

I d'c + dt,ft)| 


= g(t) 


|d e + d„ft)| 
|rd c + d„ft)| 


(B2) 


where the vertical bars indicate vector magnitude. Now define d = d c + d v and expand to 1st order in Ar 
to get 


9 fit) = g(t) 


dft)-d 
dft)-dft) 


(B3) 


Note that the term (d • d c )/(d • d) is dominated by a constant component of order d 2 /(d 2 + d 2 ) ~ 0.99, 
followed by an annually modulated term that is suppressed by a factor of order d v /d c . Thus an erroneous 
calibration dipole induces a specific error in the fit gain that can be identified and corrected for, assuming 
the time dependence of the true gain is orthogonal to this form. 
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B.2. Gain Model Fitting 

In theory, the way to do this is as follows. We have a set of data in the form of the fit gains, g/ ^ for 
each calibration sequence i, and we have a gain model, G(t;p n ), which is a function of time and a set of 
model parameters p n . Ideally we would like to fit the model to the true gain, g(t), but since we don’t know 
the true gain, the next best thing is to modify the gain model to have the same modulation form as the 
dipole gains have and to fit for this modulation simultaneously with the other gain model parameters. Thus 
X 2 takes the form 

2 [Oi - Gi(p n )] 2 [9f,i ~ Gi(p n )(l - Ar/^)] 2 

= 2- ^ ’ (B4) 

i 1 i 1 

where fd,i = (d • d c )/(d • d) evaluated at time £*, or is a function generated from simulations. 

Since the system is nonlinear, it must be minimized using a suitable nonlinear least squares routine. 
However, we can analyze the parameter covariance matrix directly by explicitly evaluating the 2nd derivative 
of x 2 with respect to the model parameters 


c -i = I d 2 X 2 

2 dpjdpk 

First compile the necessary 1st derivatives 

1 dx 2 _ [9f,i ~ Gj(Pn)( 1 ~ Arfd,j)] (Gj fd,i ) 

2 dA r cr‘ 2 

i 1 


(B5) 


(B6) 


1 d)d_ _ [gf,i - Gj(p„)( 1 - A rf dti j\ (-dGi/dp m )(l - A rf d ,i) 

2 dp m 2-* a\ 

(We evaluate the individual dG/dp m terms below.) Next the various 2nd derivatives are 

1 d 2 X 2 _ y. (Gj fd,i)(Gj fd,j) 

2 dArdAr ^ a? 

i 1 


(B7) 


(B8) 


1 d 2 X 2 

2 dArdp m 

1 g 2 X 2 = 

2 dpmdp„ ^ 


Y' (®« fd.,i)(—dGi/dpm)( 1 — Arfd,i) 
4- <T t 2 

l 1 


+ oJ2(g t - Gi), 

i 


{dGi/dp m ){\ - Arfd,i)(dGi/dp„)(l - A rf d ,i) 


+ O ^(</» — Gi). 

i 


(B9) 

(BIO) 


In the last two expressions, we neglect the term proportional to d 2 G/dp 7n dp n because the pre factor of 
(gi — Gi) is statistically zero for the least squares solution. 

Finally, we evaluate the dG/dpm terms. The gain model has the form (Jarosik et al. 2007) 


v( ti ) - vs, - /?(r RXB fe) - rg XB ) 
Tfpa{u) - r° PA 


(Bll) 


where T BXB = 290 K, and a, Vo, and Tp PA are parameters to be fit. The necessary 1st derivatives are 


dGi/da = 


V( tl )-V 0 -P(T RX B fe)-TQ XB ) 
Tfpa(^z) — ^FPA 


(B12) 
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dGi/dVo 


—a 

Tfpa(U) - T^pa ’ 


dGi/dp = 


-q(Trxb(*x) ~ Trxb) 
Tfpa(^x) ~ T f ° pa 


dGi/ dT^pk — GL 


V(tj) - Vp - 0{T RX b{U) - ^~rxb) 
(Tfpa(^i) — TOp A ) 2 


(B13) 

(B14) 
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Fig. 1. — Five-year temperature sky maps in Galactic coordinates smoothed with a 0.2° Gaussian beam, 
shown in Mollweide projection, top: K band (23 GHz), middle-left : Ka band (33 GHz), bottom-left: Q band 
(41 GHz), middle-right: V band (61 GHz), bottom-right: W band (94 GHz). 
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Fig. 2. — Five-year Stokes Q polarization sky maps in Galactic coordinates smoothed to an effective Gaussian 
beam of 2.0°, shown in Mollweide projection, top: K band (23 GHz), middle-left : Ka band (33 GHz), bottom- 
left: Q band (41 GHz), middle-right: V band (61 GHz), bottom-right: W band (94 GHz). 
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Fig. 3. — Five-year Stokes U polarization sky maps in Galactic coordinates smoothed to an effective Gaussian 
beam of 2.0°, shown in Mollweide projection, top : K band (23 GHz), middle-left : Ka band (33 GHz), bottom- 
left: Q band (41 GHz), middle-right: V band (61 GHz), bottom-right: W band (94 GHz). 
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Fig. 4. — Five-year polarization sky maps in Galactic coordinates smoothed to an effective Gaussian beam 
of 2.0°, shown in Mollweide projection. The color scale indicates polarized intensity, P = \/Q 2 + £/ 2 , and 
the line segments indicate polarization direction in pixels whose signal-to-noise exceeds 1. top: K band 
(23 GHz), middle-left : Ka band (33 GHz), bottom-left : Q band (41 GHz), middle-right : V band (61 GHz), 
bottom-right : W band (94 GHz). 
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Fig. 5. — Difference between the 5 year and 3 year temperature maps, left column : the difference in the 
maps, as delivered, save for the subtraction of a relative offset (Table 3), right column: the difference after 
correcting the 3 year maps by a scale factor that accounts for the mean gain change, ~ 0.3%, between the 
3 year and 5 year estimates, top to bottom : K, Ka, Q, V, W band. The differences before recalibration are 
dominated by galactic plane emission and a dipole residual: see Table 3, which also gives the changes for 
1 = 2,3. 
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Iteration number 



Fig. 6. — Gain convergence tests using the iterative sky map & calibration solver on a pair simulations. 
Both panels show the gain recovery as a function of iteration number for a 4-channel K band simulation. 
The initial guess was chosen to be in error by 1% to test convergence; the extrapolated output solutions 
are printed in each panel, top : Results for a noiseless simulation that includes a dipole signal (with Earth- 
velocity modulation) plus CMB and foreground anisotropy (the former is evaluated at the center frequency 
of each channel). The input gain was set to be constant in time. The extrapolated solutions agree with the 
input values to much better than 0.1%. bottom : Results for a noiseless simulation that includes only dipole 
signal (with Earth- velocity modulation) but no CMB or foreground signal. In this case the input gain was 
set up to have flight-like thermal variations. The extrapolated absolute gain recovery was in error by >0.3%, 
indicating a small residual degeneracy between the sky model and the time-dependent calibration. 
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Fig. 7. — Gain error recovery test from a flight-like simulation that includes every effect known to be 
important. Using the daily dipole gains recovered from the iterative sky map & calibration solver as input, 
the gain convergence error, shown here, is fit simultaneously with the gain model parameters, not shown, 
following the procedure outlined in Appendix B. The red trace indicates the true gain error for each WMAP 
channel, based on the known input gain and the gain solution achieved by the iterative solver on its final 
iteration. The black trace shows the gain error recovered by the fit, averaged by frequency band. The 
channel-to-channel scatter within a band is <0.1%, though the mean of Ka band error is of order 0.1%. 
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Fig. 8. — W band EE power spectrum likelihood from 1 = 2 — 7 using two separate estimation methods: 
black : maximum likelihood and red : pseudo-Cj. The vertical dashed lines indicate the best-fit model power 
spectrum based on fitting the combined Ka, Q, and V band data. The two spectrum estimates are consistent 
with each other, except at Z = 3. The maximum likelihood estimates are wider because they include cosmic 
variance whereas the pseudo-Q estimates account for noise only. Both estimates show excess power in the 
W band data relative to the best-fit model, and to the combined KaQV band spectrum, shown in Figure 6 
of Nolta et al. (2008). The extreme excess in the l = 7 pseudo-C* estimate is not so severe in the maximum 
likelihood, but both methods are still inconsistent with the best-fit model. 
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Fig. 9. — top : Simulated W band data with a large polarized thermal emission signal injected, binned by 
solar azimuth angle. The red trace shows the input waveform based on the flight mirror temperature profile 
and a model of the polarized emissivity. The black profile is the binned co-added data which follows the 
input signal very well. The thickness of the points represents the 1 a uncertainty due to white noise, bottom : 
Same as the top panel but for the 5 year flight data. The reduced \ 2 of the binned data with respect to 
zero is 2.1 for 36 degrees of freedom, but this does not account for 1/f noise, so the significance of this result 
requires further investigation. However, the much larger signal in the simulation did not produce an EE 
spectrum with features present in the flight W band EE spectrum, so the feature in the binned flight data 
cannot account for the excess / = 7 emission. 
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Fig. 10. — The EE power spectrum computed from the null sky maps, ^Ka — ^£qv, where S = Q,U are 
the polarization Stokes parameters, and S'qv is the optimal combination of the Q and V band data. The 
pink curve is the best-fit theoretical spectrum from Dunkley et al. (2008). The spectrum derived from the 
null maps is consistent with zero. 
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Optical Depth, T 

Fig. 11. — Estimates of the optical depth from a variety of data combinations. The dashed curve labeled 
Null uses the same null sky maps used in Figure 10. The optical depth obtained from Ka band data alone 
(blue) is consistent with independent estimates from the combined Q and V band data (red). The final 5 year 
analysis uses Ka, Q, and V band data combined (black). These estimates all use a 1-parameter likelihood 
estimation, holding other parameters fixed except for the fluctuation amplitude, which is adjusted to fit the 
first acoustic peak in the TT spectrum (Page et al. 2007). The degeneracy between r and other ACDM 
parameters is small: see Figure 7 of Dunkley et al. (2008). 
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Fig. 12. — The foreground-reduced Internal Linear Combination (ILC) map based on the 5 year WMAP 
data. 
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Fig. 13. — The temperature (TT) and temperature-polarization correlation (TE) power spectra based on 
the 5 year WMAP data. The addition of 2 years of data provide more sensitive measurements of the third 
peak in TT and the high-/ TE spectrum, especially the second trough. 
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Table 1. Differencing Assembly (DA) Properties 


DA 

A a 

(mm) 

i/ a 

(GHz) 

gW b 

#FWHM C 

(°) 

*o(I) d 

(mK) 

<MQ,u) d 

(mK) 

i's 6 

(GHz) 

i/ ff e 

(GHz) 

^d 6 

(GHz) 

K1 

13.17 

22.77 

1.0135 

0.807 

1.436 

1.453 

22.47 

22.52 

22.78 

Kal 

9.079 

33.02 

1.0285 

0.624 

1.470 

1.488 

32.71 

32.76 

33.02 

Qi 

7.342 

40.83 

1.0440 

0.480 

2.254 

2.278 

40.47 

40.53 

40.85 

Q2 

7.382 

40.61 

1.0435 

0.475 

2.141 

2.163 

40.27 

40.32 

40.62 

VI 

4.974 

60.27 

1.0980 

0.324 

3.314 

3.341 

59.65 

59.74 

60.29 

V2 

4.895 

61.24 

1.1010 

0.328 

2.953 

2.975 

60.60 

60.70 

61.27 

W1 

3.207 

93.49 

1.2480 

0.213 

5.899 

5.929 

92.68 

92.82 

93.59 

W2 

3.191 

93.96 

1.2505 

0.196 

6.565 

6.602 

93.34 

93.44 

94.03 

W3 

3.226 

92.92 

1.2445 

0.196 

6.926 

6.964 

92.34 

92.44 

92.98 

W4 

3.197 

93.76 

1.2495 

0.210 

6.761 

6.800 

93.04 

93.17 

93.84 


a Effective wavelength and frequency for a thermodynamic spectrum. 

b Conversion from antenna temperature to thermodynamic temperature, AT = g(v)ATA- 

c Full-width-at-half-maximum from radial profile of A- and B-side average beams. Note: beams 
are not Gaussian. 

d Noise per observation for resolution 9 and 10 /, Q, & U maps, to ~0.1% uncertainty. cr(p) = 

aJ 2 (p)- 

e Effective frequency for synchrotron (s), free- free (ff), and dust (d) emission, assuming spectral 
indices of /3 = —2.9, —2.1, +2.0, respectively, in antenna temperature units. 


Table 2. Lost and Rejected Data 


Category 

K-band 

Ka-band 

Q-band 

V-band 

W-band 

Lost or incomplete telemetry(%) 

0.12 

0.12 

0.12 

0.12 

0.12 

Spacecraft anomalies(%) 

0.44 

0.46 

0.52 

0.44 

0.48 

Planned stationkeeping maneuvers(%) 

0.39 

0.39 

0.39 

0.39 

0.39 

Planet in beam (%) 

0.11 

0.11 

0.11 

0.11 

0.11 

Total lost or rejected (%) 

1.06 

1.08 

1.14 

1.06 

1.10 
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Table 3. Change in low -l Power from 3 year Data 


Band 

1 = o a 
(mK) 

i = i a 
(mK) 

l = 2 b 
(mK 2 ) 

Z = 3 b 

(^k 2 ) 

K 

9.3 

5.1 

4.1 

0.7 

Ka 

18.9 

2.1 

2.8 

0.2 

Q 

18.3 

0.4 

2.5 

0.5 

V 

14.4 

7.3 

1.2 

0.0 

W 

16.4 

3.5 

1.0 

0.0 


a Z = 0, 1 - Amplitude in the difference map, 
outside the processing cut, in n K. 

h l — 2,3 - Power in the difference map, 
outside the processing cut, /(/ 4- 1) Cij 27r, in 
/xK 2 . 
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Table 4. Polarization \ 2 Consistency Tests a 


Multipole 

KaKa 

(10) b 

KaQ 

(50) 

KaV 

(50) 

KaW 

(100) 

QQ 

(45) 

QV 

(100) 

QW 

(200) 

vv 

(45) 

vw 

(200) 

ww 

(190) 

EE 

2 

0.727 

1.059 

1.019 

1.301 

1.586 

0.690 

1.179 

0.894 

1.078 

1.152 

3 

1.373 

0.994 

1.683 

1.355 

1.092 

1.614 

1.325 

1.005 

1.386 

1.519 

4 

1.561 

1.816 

1.341 

2.033 

0.993 

1.126 

1.581 

1.195 

1.596 

1.724 

5 

0.914 

1.313 

1.062 

1.275 

1.631 

1.052 

1.155 

0.589 

0.881 

1.252 

6 

1.003 

0.847 

0.688 

1.124 

0.740 

0.856 

1.049 

1.384 

1.168 

1.142 

7 

0.600 

0.671 

0.689 

0.936 

0.936 

0.780 

0.864 

0.900 

1.064 

1.015 

8 

1.578 

1.262 

1.337 

1.212 

1.080 

0.763 

0.608 

1.025 

0.871 

0.749 

9 

0.760 

0.710 

0.891 

0.820 

0.582 

0.726 

0.651 

0.791 

0.821 

0.795 

10 

0.494 

0.821 

0.996 

0.914 

0.656 

0.763 

0.806 

0.676 

0.891 

0.943 

EB 

2 

0.900 

1.297 

1.179 

2.074 

1.006 

0.915 

2.126 

1.242 

2.085 

2.309 

3 

0.719 

1.599 

0.651 

2.182 

1.295 

0.986 

2.739 

1.095 

3.276 

3.157 

4 

0.746 

1.702 

1.378 

1.777 

1.926 

1.110 

1.435 

1.028 

1.279 

1.861 

5 

1.161 

0.948 

0.945 

1.003 

1.149 

1.232 

1.468 

0.699 

1.122 

1.516 

6 

0.475 

1.183 

0.651 

0.687 

0.829 

1.023 

0.814 

1.201 

1.136 

0.960 

7 

1.014 

1.007 

0.829 

0.700 

0.817 

0.759 

1.112 

0.616 

0.802 

1.233 

8 

0.849 

0.897 

1.279 

0.861 

0.681 

0.689 

0.955 

1.021 

0.954 

0.996 

9 

0.743 

0.734 

1.007 

1.112 

0.820 

0.798 

0.686 

0.882 

0.808 

0.824 

10 

0.413 

1.003 

1.316 

0.859 

0.722 

0.900 

0.693 

1.124 

0.836 

0.852 

BB 

2 

2.038 

1.570 

1.244 

2.497 

1.340 

1.219 

2.529 

0.694 

1.631 

9.195 

3 

0.756 

0.868 

0.808 

1.817 

3.027 

1.717 

3.496 

0.601 

2.545 

5.997 

4 

1.058 

1.455 

1.522 

2.144 

1.007 

0.905 

1.786 

0.752 

1.403 

1.984 

5 

1.221 

1.659 

1.742 

2.036 

0.889 

1.057 

1.271 

1.078 

1.660 

1.255 

6 

0.379 

0.805 

0.483 

0.812 

1.009 

0.861 

1.238 

0.800 

0.767 

0.955 

7 

1.925 

1.594 

0.967 

1.332 

1.074 

0.817 

0.928 

0.772 

0.994 

1.024 

8 

0.804 

1.005 

0.999 

0.912 

1.069 

0.782 

0.831 

0.997 

0.879 

0.943 

9 

0.320 

0.489 

0.502 

0.450 

0.884 

0.491 

0.729 

0.748 

0.664 

0.959 

10 

1.181 

1.162 

1.028 

0.980 

1.218 

1.165 

0.951 

1.079 

0.621 

0.791 


a Table gives y 2 per degree of freedom of the independent spectrum estimates per multipole per 
band or band-pair, estimated from the template-cleaned maps. See text for details. 

b Second header row indicates the number of degrees of freedom in the reduced y 2 for that 
spectrum. See text for details. 
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Table 5. Loss Imbalance Coefficient s a 


DA 

^im,! 

(%) 

^im,2 

(%) 

K1 

0.012 

0.589 

Kal 

0.359 

0.148 

Qi 

-0.031 

0.412 

Q2 

0.691 

1.048 

VI 

0.041 

0.226 

V2 

0.404 

0.409 

W1 

0.939 

0.128 

W2 

0.601 

1.140 

W3 

-0.009 

0.497 

W4 

2.615 

1.946 


a Loss imbalance is de- 
fined as x im = ( €a — 
€b)/(€a -bee). See §6.2 
and Jarosik et al. (2007) 
for details. 
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Table 6. Cosmological Parameter Summary 


Description 

Symbol 

WMAP- only 

WMAP± BAO±SN 

Parameters for Standard ACDM Model a 


Age of universe 

to 

13.69 ±0.13 Gyr 

13.73 ±0.12 Gyr 

Hubble constant 

Ho 

71.9^2 7 km/s/Mpc 

70.1 ± 1.3 km/s/Mpc 

Baryon density 

Qb 

0.0441 ±0.0030 

0.0462 ± 0.0015 

Physical baryon density 

n b h 2 

0.02273 ± 0.00062 

0.02265 ±0.00059 

Dark matter density 

n c 

0.214 ±0.027 

0.233 ±0.013 

Physical dark matter density 

n c h 2 

0.1099 ±0.0062 

0.1143 ±0.0034 

Dark energy density 

Da 

0.742 ± 0.030 

0.721 ±0.015 

Curvature fluctuation amplitude, ko = 0.002 Mpc -1 b 

A 2 

(2.41 ±0.11) x 10“ 9 

(2.4571^) x 10- 9 

Fluctuation amplitude at Sh~ l Mpc 

O' 8 

0.796 ±0.036 

0.817 ±0.026 

i(i + i)c££/2* 

C220 

5756 ± 42 fiK 2 

5748 ± 41 fiK 2 

Scalar spectral index 

n s 

o-aealHls 

0.960i“;““ 

Redshift of matter-radiation equality 

Zeq 

3176t|U 

3280^g9 

Angular diameter distance to matter-radiation eq. c 

d>A( z e q) 

14279t}| b Mpc 

14172^^39 Mpc 

Redshift of decoupling 

Z* 

1090.51 ± 0.95 

1091.00lg;^ 

Age at decoupling 

u 

380081±|!^? yr 

375938^3^5 yr 

Angular diameter distance to decoupling c,d 

d A (z*) 

14115±i|f Mpc 

14006^ \f x Mpc 

Sound horizon at decoupling d 

r s (z *) 

146.8 ± 1.8 Mpc 

145.6 ± 1.2 Mpc 

Acoustic scale at decoupling d 

Ia{z*) 

302.08^084 

302.1ll 9 j 4 

Reionization optical depth 

T 

0.087 ±0.017 

0.084 ±0.016 

Redshift of reionization 

Z reion 

11.0 ±1.4 

10.8 ±1.4 

Age at reionization 

^reion 

427^g® Myr 

432 tlj Myr 

Parameters for Extended Models e 


Total density f 

Dtot 

1 OQQ +010 ° 
i.uyy_ 0 085 

1.0052 ± 0.0064 

Equation of state g 

W 

— 1 06 + °' 41 
i.uo_o.42 

—0 Q72+ 0 061 

u.ifi ^_ 0 060 

Tensor to scalar ratio, ko — 0.002 Mpc -1 h ' h 

r 

< 0.43 (95% CL) 

< 0.20 (95% CL) 

Running of spectral index, ko = 0.002 Mpc -1 b »* 

dn s /din k 

-0.037 ± 0.028 

— 0 032 + °' 021 

u.uoz_ 0 020 

Neutrino density •> 

PL u h 2 

< 0.014 (95% CL) 

< 0.0065 (95% CL) 

Neutrino mass j 


< 1.3 eV (95% CL) 

< 0.61 eV (95% CL) 

Number of light neutrino families k 

N eff 

> 2.3 (95% CL) 

4.4 ± 1.5 


a The parameters reported in the first section assume the 6 parameter ACDM model, first using WMAP data only 
(Dunkley et al. 2008), then using WMAP+BAO+SN data (Komatsu et al. 2008). 

h k = 0.002 Mpc -1 < * / e ff ^ 30. 

c Comoving angular diameter distance. 
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d lA(z*) = Trd A (z*)r s (z*) x . 

e The parameters reported in the second section place limits on deviations from the ACDM model, first using WMAP 
data only (Dunkley et al. 2008), then using WMAP+BAO+SN data (Komatsu et al. 2008). A complete listing of all 
parameter values and uncertainties for each of the extended models studied is available on LAMBDA. 

f Allows non-zero curvature, Qk ^ 0. 

g Allows u)/-l, but assumes w is constant. 

h Allows tensors modes but no running in scalar spectral index. 

1 Allows running in scalar spectral index but no tensor modes, 
j Allows a massive neutrino component, ^ 0. 

k Allows iV e ff number of relativistic species. The last column adds the HST prior to the other data sets. 


